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INTRODUCTION 

In  [1],  a  procedure  is  presented  for  calculating  the  electromagnetic 
field  in  the  vicinity  of  an  infinitesimally  thin  perfectly  conducting  cir¬ 
cular  cylindrical  shell  with  an  infinitely  long  slot  illuminated  by  a  TM 
wave.  A  computer  program  was  written  to  implement  this  procedure.  For  the 
shell  of  [1,  Fig.  1]  illuminated  by  the  incident  electric  field  E*  given  by 
[1.  eq.  (1)] 

Ei  _  u  ejkp  cos(<H*>  (1) 

—  — z 

the  program  calculates  the  amplitude  |e  |  of  the  z  directed  electric  field 

z 

at  equally  spaced  points  in  the  slot  aperture  and  along  the  straight  line 

that  runs  from  the  center  of  the  cylinder  to  the  center  of  the  aperture.  The 

radius  of  the  shell  is  a,  and  the  angular  width  of  the  slot  aperture  is  2<t>Q. 

The  incident  electric  field  is  the  electric  field  that  would  exist  if  the  shell 

were  absent.  In  (1),  are  cylindrical  coordinates,  and  u  is  the  unit 

— z 

vector  in  the  z-direction.  The  incident  wave  (1)  comes  from  the  direction  for 
which  <p  •  a.  In  (1),  k  »  a where  w  is  the  angular  frequency,  u  is  the 
permeability  of  the  homogeneous  medium  that  surrounds  the  shell,  and  e  is  the 
permittivity  of  this  medium. 

The  computer  program  consists  of  the  subroutines  BES,  BESJ,  BESJY, 

BESJJ,  SFEO,  DECOMP,  and  SOLVE,  and  a  main  program.  The  main  program  reads 
input  data  on  the  file  MAUTZ3.DAT  and  writes  output  data  on  file  MAUTZ6.DAT. 

The  user  who  wants  merely  to  run  the  program  is  advised  to  skip  to  Section  VIII 
and  to  read  only  the  descriptions  of  the  input  and  output  data  there. 


2 


II.  THE  SUBROUTINE  BES 

The  subroutine  BES(Y0,  N,  X,  BJ,  BY)  puts  the  Bessel  function  of  the 

first  kind  J  (X)  in  BJ(n+l)  and  the  Bessel  function  of  the  second  kind  Y  (X) 
n  n 

in  BY (n+1 )  for  n«0,l,...,N  where  X  is  a  positive  real  number.  The  subroutine 
BES  differs  from  that  in  [2,  Sec.  II]  only  in  that  common  input  variables 
Y0(1)  to  Y0(33),  PI2 ,  PI4,  PI7,  in  [2,  p.  6]  have  been  changed  to  the  sub¬ 
routine  arguments  Y0(1)  to  Y0(33).  Y0(34) ,  Y0(35) ,  and  Y0(36),  respectively, 
and  the  common  output  Variables  BJ  and  BY  have  been  changed  to  the  subroutine 
arguments  BJ  and  BY,  respectively.  The  symbol  0  in  Y0  denotes  zero  rather 
than  the  letter  0.  In  other  instances  when  it  is  obvious  from  the  context 
that  a  zero  is  meant,  the  slash  is  not  used. 

Minimum  allocations  are  given  by 

DIMENSION  BJ  (N+1) ,  BY(N+1),  AJ(114-N+2*[X] )  where  [X]  is  the  largest 


integer  that  does  not  exceed  X. 
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LISTING  OF  THE  SUBROUTINE  BES 
SUBROUTINE  BESlYO.N.X.BJ.BY) 

D I  HENS I ON  Y0136) .BJ 1 100) .BY ( 100) . AJ 1 100) 
HZ*10+N+2«1FIX1X) 

IF(X.LT. l.E-3)  HZ-4+N 
AJ1HZ+1)*0. 

AJlHZ)*l.E-20 

Hl-HZ-1 

X2-2./X 

DO  16  1*1, Hi 

HK-HZ-K 

AJ  <  HI ) «X2»FL0AT (HI) »AJ 1HK+1 ) -A J ( HK+2 > 

16  CONTINUE 
ALP-.5»AJ11) 

DO  15  J*3.HZ.2  53  Sl*Y012 

ALP-ALP+AJ ( J )  54  S3«Z» < 1 

15  CONTINUE  55  S5-X-YC 

ALP *2. *ALP  55  BT( 1 )*S 

JJP*11*!  57  IF1N.LI 

DO  17  1*1,  HP  55  S2-T0C 

BJ(K)*AJ(K)/ALP  59  S4-Z*n 

17  CONTINUE  60  S8-X-TC 

IF<X-8. >18.18. 19  51  BY12)*S 


S 1  -YO 1 22 )  ♦  YO 1 23 )  «Z  1  ♦  YO 1 24 ) «  Z2 
S3*Z* 1 YO 1 28 ) +Y0 1 29 ) *Z 1 + YO 1 30 ) *Z2 ) 
S5-X-Y0135) 

BY ( 1  )*S»  <S1*SIN(S5) +S3»C0S(S5) > 

IF(N.LE.O)  RETURN 

S2-Y0  1 25  > ♦ YO 1 26 ) »Z1 ♦ YO 1 27  >  »Z2 

S4*Z» 1Y0131 )+Y0t32>  »Z1+Y0133) »Z2) 

S6*X-Y0136) 

BY12)*S* 1S2*SIN1S6) +S4«C0S1S6) ) 


18  Z1-X*X  62  20  IF1NP.LE.2)  RETURN 

5^0*  63  DO  21  1*3.  HP 

J3*0*  64  K2-K-2 

85  BY 1 I) *X2*FL0AT  1 12 ) *BY 1 I- 1 ) -BY 1 12 ) 

“*°*  66  21  CONTINUE 

Z5*°*  67  RETURN 

27-0  •  68  END 

IF1Z1.LT. l.E-11)  GO  TO  22 
Z2*Z1»Z1 
Z3*Z2«Z1 

IF1Z1.LT. l.E-6)  GO  TO  22 
Z4*Z3«Z1 
Z5*Z4»Z1 

IF1Z1.LT. l.E-4)  GO  TO  22 
Z6*Z5«Z1 
Z7*Z8»Z1 

22  S-l Y01 1 >+Y012> »Z1+Y013) «Z2+Y0<4> *Z3*Y0<5>  tZ4+Y016) *Z5+ 

1Y017) »Z8+Y018) »Z7)/ 1 Y019)*Y01 10) tZl+YOl 1 1 ) »Z2) 

TLG-ALOG(X) 

BY 1 l ) *S+YO 134 ) «BJ 1 1 ) #TLG 
IF1N.LE.0)  RETURN 

S*(TO( 12) ♦YOl 13) »Z1+Y0( 14) »Z2*Y01 15) *Z3*Y0< 16) *Z4* 

1Y01 17) »Z5+Y01 18) »Z6) / 1Y01 19) +Y0120) »Z1+Y0<21 ) *Z2) 

BY 1 2 ) *X»S+ YO 1 34 ) *  1 B J 1 2 ) »TLG- 1 . / X ) 

GO  TO  20 

19  Z*8./X 
Z1-Z»Z 
Z2*Z1*Z1 

S«SflRT 1 TO 1 34 ) / X ) 


P& 


III.  THE  SUBROUTINE  BESJ 

The  subroutine  BESJ(X,  BJZ,  BJ1)  puts  Jq(X)  in  BJZ  and  J^(X)  in 

BJ1  where  X  is  a  positive  real  number. 

If  X  <  3,  then  J  (X)  and  J. (X)  are  calculated  in  lines  5  to  20  in 
o  l 

the  listing  of  BESJ  at  the  end  of  this  section.  As  suggested  in  [3,  Sec. 

9.12.,  Example  1],  J^(X)  and  are  set  equal  to  the  arbitrary 

-20 

values  of  zero  and  10  ,  respectively,  where 

MZ  -  11  +  2[X]  (2) 

where  [X]  is  the  largest  integer  that  does  not  exceed  X.  This  is  done  in 
lines  5  to  7.  In  DO  loop  16,  the  recurrence  relation  [3,  Eq.  (9.1.27)] 

Vx>  -  Wx>  '  J„+2(x)  «> 

is  used  to  calculate  (J  (X),  n  »  MZ-2,  MZ-3,.  ..,0}.  Line  12  puts  J__  ,  (X) 

n  MK— 1 

in  BJ(MK).  According  to  [3,  eq.  (9.1.46)],  the  calculated  values  of  J  (X) 

o 

and  J^(X)  have  to  be  normalized  by  dividing  by  a  where 

a  -  J  (X)  +  2J,(X)  +  2J. (X)  +  ...  (4) 

The  division  by  a  is  performed  in  lines  19  and  20. 

If  X  >  3,  then  J  (X)  and  J, (X)  are  calculated  in  lines  22  to  38. 

o  1 

The  polynomial  approximations  [3,  Eqs.  9.4.3  and  9.4.6]  are  used.  With 
regard  to  [3,  Eq.  9.4.3],  lines  28  and  29  put  f  in  FZ,  lines  30  and  31 

put  9q  in  TZ,  and  line  37  puts  Jq(X)  in  BJZ.  With  regard  to  [3,  Eq.  9.4.6], 

lines  32  and  33  put  f^  in  FI,  lines  34  and  35  put  0^  in  Tl,  and  line  38 


1C  LISTING  OF  THE  SUBROUTINE  BESJ 

2  SUBROUTINE  BESJ(X,BJZ,BJ1) 

3  01 HENS I ON  BJ ( 18) 

4  IF(X-3. )  17,17,18 

5  17  HZ-11+2MFIXIX) 

6  BJ<HZ+l>-0. 

7  BJ(HZ)«l.E-20 

8  Hl-HZ-1 

9  X2-2./X 

10  DO  18  K«1,H1 

11  HK-HZ-K 

12  B J ( HE ) «X2»FL0 AT ( HK ) *B J  <  HK+ 1 ) -B J ( HK+2 ) 

13  16  CONTINUE 

14  ALP«.5«BJ<1> 

15  DO  19  J-3.HZ.2 

16  ALP-ALP+BJ(J) 

17  19  CONTINUE 

18  ALP-2. »ALP 

19  BJZ«BJ( 1)/ALP 

20  BJl»BJ<2)/ALP 

21  RETURN 

22  18  X1-3./X 

23  X2-X1«X1 

24  X3«X2*X1 

25  X4-X3«X1 

26  X5-X4*X1 

27  X6-X5»X1 

28  FZ- . 7978846- . 77E-6 »  X 1  - . 55274E-2*  X2- . 95 1 2E-4  *  X3+ 

29  1 . 137237E-2*X4~. 72805E-3»X5+. 14476E-3»X6 

30  TZ-X-. 7853982-. 4166397E-1 »X1-. 3954E-4»X2+. 262573E-2»X3- 

31  1 . 54125E-3»X4-. 29333E-3»X5+. 13558E-3«X6 

32  F1-. 7978846+. 156E-05«X1+. 1659667E-01«X2+. 17105E-03*X3- 

33  1 . 24951 1E-02»X4+. 1 13653E-02»X5-. 20033E-03*X6 

34  T 1  «X~2 . 356 1 94+ .  1 24996 1  »X  1  ♦ .  565E-04  •  X2- .  637879E-02*  X3+ 

35  l. 74348E-03»X4+. 79824E-03«X5-. 29166E-03*X6 

36  XX-SQRT(X) 

37  BJZ-FZ/XX«COS<TZ> 

38  BJ1»F1/XX»C0S<T1) 

39  RETURN 

40  END 
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IV.  THE  SUBROUTINE  BESJY 


The  subroutine  BESJY(N,  X,  BJY)  puts  J„(X)  Y„(X)  in  BJY  where 
JN  and  Yn  are,  respectively,  the  Bessel  functions  of  the  first  and 
second  kinds  of  order  N.  Here,  X  is  a  non-negative  real  number,  and 
N  is  large  enough  so  that  Debye's  asymptotic  expansions  [3,  eqs.  (9.3.7) 
and  (9.3.8)]  apply  to  J^(X)  and  YN(X).  The  subroutine  BESJY  is  that  in 
[2,  Sec.  VI]  with  DO  loop  11  removed,  the  arguments  N1  and  N2  replaced 


by  N,  and  the  array  BJY  replaced  by  the  single  variable  BJY. 


LISTING  OF  THE  SUBROUTINE  BESJY 
SUBROUTINE  BESJY (N.X. BJY) 

F1«N 

F2-F1»F1 

F3«F2»F1 

F4«F3»F1 

XN-X/F1 

T1*1./SQRT(1.-XN«XN) 

T2«T1»T1 

T3«T2»T1 

T4«T3»Ti 

T5«T4»T1 

TB«T5*T1 

T7-T6»T1 

T8«T7»T1 

T9»T8»TI 

T10»T9«T1 

T12«T10»T2 

U1«(3.»TI-5.»T3)/(24.«F1) 

U2»(81. *T2-462.»T4+385. *T6)/ ( 1152. »F2) 

U3» (30375. »T3- 389603. *T5*765765. *T7-425425. *T9)/ <414720. *F3) 
U4- ( 4485 1 25. »T4-94 1 2 1 68 . E+0 1 *T6*3499224 . E*02*T8-446 1 857 . E+02« 
1T10+ 1859 107. E+02*T12)/ (3981312. E*0i*F4) 

U5-1.+U2+U4 

U6-U1+U3 

BJY— Tl  /  (3. 141593»F1 ) »  (U5+U6)  *  (U5-U6) 

RETURN 

END 


V.  THE  SUBROUTINE  BESJJ 


The  subroutine  BESJJ (N,  NR,  X, 


Jn(X(L)) 

BJJ)  puts  -j — (xTNRT)  ln  BJJ^  for 


L  *  1,2 . NR.  Here,  X(L)  _>  0  and  N  is  large  enough  so  that  Debye's 

asymptotic  expansion  applies  to  J^.  Minimum  allocations  are  given  by 

DIMENSION  X(NR) ,  BJJ(NR),  U5(NR),  U6(NR) 

Consider  JN(XX)  where 


XX  -  X(L)  (5) 

Debye's  asymptotic  expansion  for  J.,(XX)  is  [3,  eq.  (9.3.7)] 

N 

N(tanha-a)  4  u  (coth  a) 

JN(XX)  -  - - (1  +  I  r - )  (6) 

’  /2trN  tanh  a  k**l  N 


where  u^  is  given  by  [3,  formula  9.3.9]  and 


.  r.  io ri 

tanh  a  -  yl  -  (— ) 


(7) 


From  [4,  formula  702.],  we  have 


1  .  ,  1  +  tanh  a. 

a  *  —  ln  (- - - — ) 

2  1  -  tanh  a 

Substituting  (7)  into  (8)  and  rationalizing  the  denominator  of  the 
argument  of  the  logarithm,  we  obtain 

N 

a  *  ln  [—  (1  +  tanh  a)] 

Substitution  of  (9)  into  (6)  gives 

tanh  a 


J  (XX)  *  _ - _  S. 

N  /5Vn  tanh""a  N<1+tilnh 


7  ,S  (1+  .I 


4  u^(coth  a) 


(8) 


(9) 


(10) 


k=l  N 


Thanks  to  (10),  we  obtain 


JN(XX)  =  U6(L)  U3(L)  N 

J  (X(NR))  =  'U6 (NR)  U5 (NR) 


(ID 
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where 


U5(L) 


XX  e 


tanh  a 


1  +  tanh  a 


(12) 


and 

4  u,  (coth  a) 

1+  I  - 

k»l  N 

U6(L)  -  - -  -  ---- — -  (13) 

/tanh  a 

Here,  tanh  a  depends  on  L  by  means  of  (7)  with  XX  given  by  (5). 

In  the  listing  of  the  subroutine  BESJJ  at  the  end  of  this  section, 
line  5  inside  DO  loop  11  obtains  XX  of  (5).  Lines  23  to  26  put  u^,  u^,  u^» 
and  u^  of  (13)  in  Ul,  U2,  U3,  and  U4,  respectively.  Line  28  obtains  U5(L) 
of  (12).  Line  29  obtains  U6(L)  of  (13).  Inside  DO  loop  12,  line  34  puts 
Jn(XX)/Jn(X(NR))  of  (11)  in  BJJ(L). 


.  ** 


1C  LISTING  OF  THE  SUBROUTINE  BESJJ 

2  SUBROUTINE  BESJJ <N. NR. X, BJJ) 

3  DIMENSION  X< 1 1) .BJJ ( 11) ,U5< 11) .UB< 11 > 

4  00  11  L-l.NR 

5  XX-X(L) 

6  Fl-N 

7  F2»F1#F1 

8  F3*F2»F1 

9  F4«F3»F1 

10  XN-XX/F1 

11  TT«SQRT(1.-XN«XN) 

12  Tl-l./TT 

13  T2*T1»T1 

14  T3»T2»T1 

15  T4-T3»T1 

16  T5*T4«T1 

17  T6«T5*T1 

18  T7«T8*T1 

19  T8»T7»T1 

20  T9»T8»T1 

21  T10«T9»T1 

22  T12«T10»T2 

23  U1«<3.»T1-5.»T3)/(24.»F1) 

24  U2a(81« •T2-482. »T4+385. »T8)/ ( 1 152. »F2> 

25  U3«< 30375. »T3-369603. •T5+765785. »T7-425425. *T9)/ (414720. »F3> 

26  04^(4465125. »T4-9412l68.E+01»T6+3499224. E+02»T8-4461857. E+02« 

27  1T10+1859107.E+02«T12)/ <3981312. E+0l«F4) 

28  U5(L)*XX*EXP(TT)/ < 1. ♦TT) 

29  U6<L)»( 1. ♦U1+U2+U3+U4)/SQRT(TT) 

30  11  CONTINUE 

31  S5«U5(NR) 

32  S6«U6<NR> 

33  DO  12  L-t.NR 

34  BJJ (L)sU6(L)/S6* (U5(L)/S5) **N 

35  12  CONTINUE 

36  BJJ(NR)>1. 

37  RETURN 

38  END 
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VI.  THE  SUBROUTINE  SFEO 

The  subroutine  SFEO(N,  P,  SE»  SO,  FE,  FO)  puts  of  [1,  eqs. 
(A-18)  to  (A-21) ]  or  [1,  eqs.  (A-32)  to  (A-35) ]  in  SE(J),  of  [1,  eqs. 
(A-22)  to  (A-25) ]  or  [1,  eqs.  (A-37)  to  (A-40)j  in  S0(J),  F^  of  [1,  eqs. 
(B-5)  to  (B-8) ]  or  [1,  eqs.  (B-18)  to  (B-21)]  in  FE(J),  and  F^  of  [1, 
eqs.  (B-9)  to  (B-12)]  or  [1,  eqs.  (B-22)  to  (B-25)]  in  F0(J),  all  for 
J-1,2,3,  and  4.  The  variables  N  and  P  are  input  arguments  defined  by 


where  n  and <Pq  appear  in  [1,  eq.  (B-4)]  so  that  [1,  eq.  (B-4)]  becomes 

b  -  N*P  (15) 

The  subroutine  SFEO  calls  the  subroutine  BESJ. 

The  subroutine  SFEO  is  listed  at  the  end  of  this  section.  If 
b  >  2,  line  11  of  this  listing  puts  JQ(b)  and  J^(b)  in  BJZ  and  BJ1, 
respectively.  Lines  15  to  19  put  S^  of  [1,  eqs.  (A-18)  to  (A-21)]  in 
SE(J)  for  J-1,2,3,  and  4.  Lines  22  to  26  put  S^  of  [1,  eqs.  (A-22)  to 
(A-25)]  in  SO(J)  for  J-1,2,3,  and  4.  Lines  30  to  33  put  F^  of  [1,  eqs. 
(B-5)  to  (B-8)]  in  FE(J)  for  J-1,2,3,  and  4.  Lines  36  to  40  put  F^  of 
[1,  eqs.  (B-9)  to  (B-12)]  in  F0(J)  for  J-1,2,3,  and  4. 

If  b  <  2,  lines  42  and  43  put  e  <J>  / 2  in  PP.  Lines  46  to  53  put 
S^  of  [1,  eqs.  (A-32)  to  (A-35)]  in  SE(J)  for  J-1,2,3,  and  4.  Lines  55 
to  62  put  S^  of  [1,  eqs.  (A-37)  to  (A-40)]  in  S0(J)  for  J-1,2,3,  and  4. 
Lines  63  to  70  put  F^  of  [1,  eqs.  (B-18)  to  (B-21)]  in  FE(J)  for  J-1,2,3, 
and  4.  Lines  71  to  78  put  F^  of  [1,  eqs.  (B-22)  to  (B-25)]  in  F0(J)  for 
J-1,2,3,  and  4. 


1C  LISTING  OF  THE  SUBROUTINE  SFEO 

2C  THE  SUBROUTINE  SFEO  CALLS  THE  SUBROUTINE  BESJ. 

3  SUBROUTINE  SFEOCN.P.SE.SO.FE.FO) 

4  D I HENS I ON  SE(4> ,S0<4> .FE<4) . F0(4) 

5  FN-N 

6  B»FN*P 

7  B2aB*B 

8  B4»B2*B2 

9  B6»B4»B2 

10  IFIB.LE.2. )  GO  TO  11 

11  CALL  BESJCB.BJZ.BJl) 

12  SZ-BJZ/FN 

13  BZaB«SZ 

14  B1-BJ1/FN 

15  SE( 1 )aBl 

16  SE(2)«( (B2-6. >  *Bl+3. *BZ)/B2 

17  SE(3)a< (B4-27. *B2+120. ) »Bl+6. * (B2-10. ) *BZ)/B4 

18  SE<4>«< (B6-63. *B4+1200. *B2-5040. ) *Bl+3. • (3. *B4-95. »B2+840.  ) *B 

19  1Z)/B6 

20  BZ»SZ 

21  Bl-Bl/B 

22  SOd)"^.  *B1-BZ 

23  S0(2)a< (5.*B2-24. )»B1-(B2-12. )*BZ)/B2 

24  SO  <3>a< (8. *B4-168. »B2+720. ) »Bl-(B4-39. *B2+360. ) *BZ>/B4 

25  SO (4)-( ( 1 1 . «B6-537. *84+9720. *82-40320. ) *B1- (B6-81 . »B4+2340. *B 

26  12-20160. )»BZ)/B6 

27  SN-SIN<B)/FN 

26  SSaCOS(B)/FN 

29  CS>B«SS 

30  FEdl-SN 

31  FE(2)a<2. »CS+<B2-2. ) *SN)/B2 

32  FE(3)»(4.*<B2-6. >*CS+<B4-12. «B2+24. )«SN)/B4 

33  FE(4)>(6. * (B4-20. *B2+120. ) »CS+ (B6-30. »B4+360. *B2-720.  ) *SN)/B6 

34  SN-SN/B 

35  CSaSS 

36  FO ( 1 ) -SN-CS 

37  FO ( 2 ) a  <  3. » ( B2-2. ) *SN- <  B2-6. >  *CS ) / B2 

38  F0<3>a<5.»(B4-12.*B2+24. > *SN-(B4-20. *B2+120. >*CS>/B4 

39  F0(4)«(7. » (B6-30. «B4+380. »B2-720. )»SN-(B6-42. »B4+840. »B2-5040 

40  l.)»CS)/B6 

41  RETURN 

42  1 1  PP-P 

43  IF(N.EQ.O)  PP-.5*PP 

44  68*86*82 

45  B10aB8»B2 

46  SE( 1 )«PP/2. * ( 1. -B2/8. +84/ 192. -B6/9216. ♦BS/ 737280. -BIO/ 884736. 

47  1E+02) 

48  SE(2)-PP/8. * ( 1 . -B2/4. +B4*5. /384. -BB*7./ 23040. *88/245760. -BIO* 

49  1 11./ 3096576. E+02) 

50  SE(3) »PP/ 16. *<1. -B2*5. / 16. +B4*7. /384. -B6*7. / 15360. +B6* 1 1 . / 1 72 

51  l 0320. -BIO* 143./ 247726 l.E+03) 

52  SE (4) «PP*5. / 128. * ( 1 . -B2»7. /20. +B4*7. /320. -B6* 11./ 19200. +B8* 14 

53  l 3. / 1 72032 . E+02-B 1 0* 1 43 . / 1 857946 . E+03  > 


/. 
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54  PBaP»B 

55  SO  ( 1  )»PB/8.  *  ( 1.  -B2/ 12. +B4/384.  -B6/23040.  +B8/221 1840.  -B10/3096 

56  1576. E+02) 

57  SK)(2)aPB/16. » ( l.-B2»5./48. +B4»7./ 1920. -B6/ 15360. +B6«1 1./ 15462 

58  188. E+01-B10» 13. /2477261 . E+03 ) 

59  SO (3) >PB«5. / 128. Ml. -B2»7. /60. +B4»7. / 1600. -B6« 11./ 134400. +B6» 

60  1 143. / 1548288. E+02-B10* 13. / 1857946. E+03) 

61  SO  <4)  «PB»7. /256. « ( 1 . -B2/8. +B4« 1 1 . /2240. -B6* 143/ 1505280. +B8« 14 

62  13. / 1300562. E+02-B 10» 17 . / 2000864 . E+03 ) 

63  FE<  1  >  aPM  1 . -B2/6. +B4/ 120. -B6/5040. +B6/362860. -BIO/ 399 168. E+02 

64  1) 

65  FE ( 2 ) «P/ 3. » < 1 . -B2* . 3+B4/ 56 . -B6/ 2 1 60 . +B8/ 1 47640. -B 1 0/ 1 57248 . E+ 

66  102) 

67  FE(3)«P/5.M1.-B2»5./ 14. +B4»5./216. -B6/ 1584. +B8/ 104832. -BIO/ 1 

68  108864. E+02) 

69  FE(4>aP/7. Ml. -B2«7. / 18.+B4»7. /264. -B6«7. /9360. +B6/86400. -BIO 

70  1/8812800. ) 

71  FQ< l)aPB/3. M 1.-B2*. 1+B4/280.-B6/ 15120. +B6/1 330560 . -B 1 0/ 1 7297 

72  128. E+02) 

73  FO  <2)«PB/5. Ml. -B2»5./42. +B4/216. -B6/ 1 1088. +B8/943488. -BIO/ 1 1 

74  197504. E+02) 

75  FO (3)«PB/7. • < 1 . -B2«7. /54. +B4»7. / 1320. -B6/9360. +B8/777600. -BIO 

76  1/969408. E+02) 

77  F0(4)aPB/9. « ( 1 . -B2*3./22. +B4»3. /520. -B6/6400. +B8/685440. -BIO/ 

78  1842686. E+02) 

79  RETURN 

80  END 


VII.  THE  SUBROUTINES  DECOMP  AND  SOLVE 

The  subroutines  DECOMP(N,  IPS,  UL)  and  SOLVE(N,  IPS,  UL,  B,  X) 
solve  a  system  of  N  linear  equations  in  N  unknowns.  The  input  to 
DECOMP  consists  of  N  and  the  N  by  N  matrix  of  coefficients  on  the 
left-hand  side  of  the  matrix  equation  stored  by  columns  in  UL.  The 
output  from  DECOMP  is  IPS  and  UL.  This  output  is  fed  into  solve.  The 
rest  of  the  input  to  solve  consists  of  N  and  the  column  of  coefficients 
on  the  right-hand  side  of  the  matrix  equation  stored  in  B.  SOLVE  puts 
the  solution  to  the  matrix  equation  in  X. 

Minimum  allocations  are  given  by 

COMPLEX  UL(N*N) 

DIMENSION  SCL(N) ,  IPS(N) 
in  DECOMP  and  by 

COMPLEX  UL(N*N),  B(N) ,  X(N) 

DIMENSION  IPS(N) 

in  SOLVE 

More  detail  concerning  DECOMP  and  SOLVE  is  on  pages  46-49  of  [5]. 
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LISTING  OF  THE  SUBROUTINE  DECOHF 
SUBROUTINE  DECOHPIN. IPS.UL) 

CO HP LEX  UL( 16). PIVOT. EH 
D I HENS I ON  SCL(4), IPS<4> 

00  5  1-1. N 
IPS< I )-l 
RN«0. 

Jl-I 

DO  2  J-l.N 

ULH-ABS ( REAL ( UL ( J 1 > ) ) +ABS ( A I  HAG ( UL ( J 1 ) > ) 
Jl-Jl+N 

IF(RN-ULH)  1.2.2 

RN-ULH 

CONTINUE 

SCL(  D-l./RM 

CONTINUE 

NH1-N-1 

K2-0 


00  17  E-1.NH1 
BIG-O. 

DO  11  I-K.N 
IP-IPS(I) 

IPE-IP+K2 

S I ZE- ( ABS ( REAL ( UL ( IPK) ) ) ♦ABS(AIHAG(UL( IPK) ) > ) »SCL< IP) 
IF(SIZE-BIG)  11.11.10 


BIG-SIZE 

IPV-I 

CONTINUE 

IF(IPV-K)  14.15.14 
J-IPS(K) 

IPS(K)-IPSdPV) 

IPS(IPV)-J 
KPP- 1  PS  ( K )  *K2 
P I VOT-UL ( KPP ) 

KP1-K+1 
DO  16  l-KPl.N 
EP-KPP 
IP-IPS( I ) +K2 
EH— UL(IP)/PIVOT 
UL  ( I P )  —EH 
00  16  J-KP1.N 
IP-IP+N 
KP-KP+N 

UL( I P ) -UL ( IP)+EH*UL<KP) 

CONTINUE 

K2-K2+N 

CONTINUE 

RETURN 

END 

LISTING  OF  THE  SUBROUTINE  SOLVE 
SUBROUTINE  SOLVEIN. IPS.UL.B.X) 
CO HP LEX  UL(16).B(4).X(4).SUH 
0 1 HENS I ON  1  PS ( 4 ) 


54 

NP1-N+1 

55 

IP-IPS(l) 

56 

X ( 1 > -B ( IP) 

57 

00  2  1-2. N 

56 

IP-IPS(I) 

59 

IPB-1P 

60 

IHi-l-1 

61 

SUH-O. 

62 

00  1  J-1.IH1 

63 

SUH-SUH+UH  IP)  *X(  J> 

64 

1 

IP-IP+N 

65 

2 

X ( J  > -B ( 1 PB  > -SUH 

66 

K2-N*(N-1> 

67 

IP-IPS(N) ♦K2 

66 

X(N)-X(N)/UL( IP) 

69 

00  4  I  BACK-2. N 

70 

I -NP 1 - 1  BACK 

71 

K2-K2-N 

72 

IPt«IPS< I ) ♦K2 

73 

IP1-IM 

74 

SUH-O. 

75 

IP-IPI 

76 

00  3  J-IP1.N 

77 

IP-IP+N 

78 

3 

SUH-SUH^UL ( IP)-X(J) 

79 

4 

X( I )-(X< I )-SUH)/UL( 

60 

RETURN 

61 

END 

1 


& 
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VIII.  THE  MAIN  PROGRAM 

The  main  program  accepts  input  data  from  the  file  MAUTZ3.DAT  and 
writes  output  data  on  the  file  MAUTZ6.DAT.  See  the  open  statements  on 
lines  10  and  11  of  the  listing  of  the  main  program  at  the  end  of  this  sec¬ 
tion.  The  main  program  calls  the  subroutines  BES,  BESJY,  BESJJ,  SFEO, 
DECOMP,  and  SOLVE.  The  subroutine  SFEO  calls  the  subroutine  BESJ.  In 
this  section,  first  the  input  data  are  described,  then  the  output  data 
are  described,  relevant  formulas  are  given,  the  main  program  is  de¬ 
scribed  verbally,  and  finally  the  main  program  and  sample  input  and  output 
data  are  listed. 

The  input  data  are  read  early  in  the  main  program  from  the  file 
MAUTZ3.DAT  according  to 

READ(2O,25)(Y0(I),  I  -  1,  36) 

25  FORMAT (5E14. 7) 

READ(20 ,28)  N1 ,  N2 ,  N3,  NA,  NR 
28  FORMAT (13,  15,  313) 

DO  31  JW  -  1,  N3 

READ(20,32)  NN,  IP,  IB,  X,  P,  ALP 
32  FORMAT (3 13,  3E14.7) 

31  CONTINUE 

The  input  array  Y0  should  not  deeply  concern  the  user  because  he  will  never 
have  to  change  it.  The  values  of  the  elements  of  Y 0  appear  in  the  sample 
data  listed  after  the  main  program.  The  curious  reader  will  find  the 
values  of  Y0(1)  to  Y0(33)  explained  on  pages  4  and  5  of  [2].  The  values 
of  Y0(34),  Y0(35),  and  Y0(36)  are  2/",  n/4,  and  3',/4,  respectively. 


^  A -*•  *»*-**'*-'•  -*’  -Vlvl'y  >>  -V  -**  -**  -V-V. 


The  Bessel  functions  in  the  summations  in  [1,  eqs.  (41),  (42), 

(47),  and  (48)]  are  approximated  by  Debye's  asymptotic  expansions 
whenever  N1  <_  n  £  N2,  and  these  summations  are  truncated  at  n  ■  N2. 

The  summations  in  [1,  eqs.  (44)  and  (50)]  are  truncated  at  n  *  Nl-1. 

For  the  z  directed  electric  field  E^(-M)  along  the  straight  line  that 
runs  from  the  center  of  the. cylinder  to  the  center  of  the  aperture,  the 
Bessel  functions  in  the  first  summation  in  (62)  are  approximated  by 
Debye's  asymptotic  expansions  whenever  N1  <_  n  <_  N2,  and  this  summation 
is  truncated  at  n  *  N2. 

The  aperture  field  is  evaluated  at  NA  equally  spaced  points  in 

the  aperture.  .iere,  NA  >_  2.  With  reference  to  [1,  Fig.  1],  the  first 

point  is  at  the  beginning  of  the  aperture  at  $  ■  an<*  t*ie  point 

is  at  the  end  of  the  aperture  at  <J>  *  4>  •  The  interior  field  Eb(-M)  is 

o  z  — 

evaluated  at  NR  equally  spaced  points  along  the  line  that  runs  from  the 
center  of  the  cylinder  to  the  center  of  the  aperture.  Here,  NR  _>  2. 

The  first  point  is  at  the  center  of  the  cylinder,  and  the  NRth  point  is 
at  the  center  of  the  aperture. 

The  NA  evaluations  of  the  aperture  field  and  the  NR  evaluations 
of  the  interior  field  Eb(-M)  are  done  inside  DO  loop  31,  that  is,  these 
evaluations  are  done  N3  times,  once  for  each  value  of  JW.  If  NN  *  4, 
these  fields  are  obtained  by  solving  [1,  eqs.  (40)  and  (46)]  in  which  Ye 
and  Y°  are,  as  stated  in  [1],  4x^  matrices.  If  1  NN  <_  3,  the  fields 
are  obtained  by  solving  [1,  eqs.  (40)  and  (46)]  with  Y  replaced  by  the 
NN x NN  matrix  in  its  upper  left-hand  corner,  with  Y°  replaced  by  the 
NN x NN  matrix  in  its  upper  left-hand  corner,  and  with  Ve ,  T6 ,  and  V° ,  and  I 
truncated  accordingly.  We  assume  that  1  NN  _<  4. 


In  DO  loop  31,  IP  *  p  where  p  is  the  non-negative  integer  upon 
which  the  alternative  expansion  functions  [1,  eqs.  (30)  and  (34)]  de¬ 
pend.  The  main  program  was  written  under  the  assumption  that  IP  <  Nl. 


If  there  is  an  integer  p  such  that  |  ( k.a )  |  is  very  small,  then  the  user 

should  take  IP  *  p.  Here,  k  is  the  wave  number  and  a  is  the  radius  of 

the  cylinder.  However,  if  there  is  no  value  of  p  such  that  |J  (ka) I  is 

P 

very  small,  a  non-negative  value  of  IP  must  still  be  chosen,  and  the 
alternative  expansion  functions  will  be  given  by  [1,  eqs.  (30)  and  (34)] 
with  p  ■  IP.  In  this  case  the  fields  calculated  by  the  procedure  of  [1] 
should  not  depend  on  p  so  that  IP  can  be  any  non-negative  integer  less 
than  Nl.  If  IB  ^  0,  then  the  main  program  uses  the  value  of  Jpp(ka) 
calculated  by  the  subroutine  BES.  If  IB  *  0,  then  Jpp(ka)  is  set  equal 
to  zero.  If  the  circular  cylinder  of  radius  a  was  resonant  such  that 
Jjp(ka)  *  0,  then,  without  IB,  it  would  be  virtually  impossible  to  obtain 
Jpp(ka)  *  0  because  the  exact  resonant  value  of  ka  would  be  difficult  to 
enter,  and  even  if  it  could  be  entered,  the  subroutine  BES  would  probably 
not  give  exactly  zero  for  Jpp(ka). 

In  DO  loop  31, 


where  p,  ka,  and  a  appear  in  [1].  Here,  P  and  ALP  are  in  radians. 

The  meaning  of  p  was  clarified  in  the  previous  paragraph,  ka  is  the  prod¬ 
uct  of  the  wave  number  k  and  the  radius  a  of  the  cylinder,  <t  is  one  half 

o 

the  angular  width  of  the  aperture,  and  *  is  the  value  of  i  in  the  direc- 
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tion  from  which  the  incident  wave  comes.  The  angles  <t>Q  and  a  are  shown 
in  [1,  Fig.  1]. 

Minimum  allocations  are  given  by 

DIMENSION  BJ (MAX (N1*NR ,  NA)),  BY (N1*NR) ,  XR(NR) ,  RP(NR) 

DIMENSION  RN(NN*NR) ,  RE (NR) ,  BJJ(NR) 

where  MAX  denotes  the  maximum  value. 

The  main  program  writes  output  data  on  the  file  MAUTZ6.DAT.  The 
input  data  of  the  first  two  read  statements  are  written  out  immediately 
after  they  are  read  in. 

The  following  write  statements  are  in  DO  loop  31: 

WRITE (21 ,  27)  NN,  IP,  IB,  X,  P,  ALP 
27  FORMAT ( '  NN-' ,  13,',  IP»’ ,  13,  ',  IB  -  ',  13/ 

1*  X  -  ’,  E14.7,  *,  P-',  E14.7,  ',  ALP-*,  E14.7) 

WRITE(21 ,  33)  (BJ(NP),  NP-JJ,  JJN) 

33  FORMAT ( '  BJ'/(1X,  5E14.7)) 

WRITE(21 ,  19)  (BY(NP),  NP  =  JJ,  JJN) 

19  FORMAT('  BY'/(1X,  5E14.7)) 

WRITE(21 ,  24)  (SE(J) ,  J«1 ,  NN) 

WRITE(21 ,  24)  (SO(J),  J«1 ,  NN) 

24  FORMAT (IX,  4E14.7) 

WRITE(21,  49)  U2 ,  SI 

49  FORMAT ( '  E-' ,  2E14.7,',  ABS(E)  *  E14.7) 

WRITE(21 ,64)  (BJ (L) ,  L  »  1,  NA) 

64  FORMAT ('  APERTURE  FIELD  AMPLITUDE' /(IX,  4E14.7)) 

WRITE(21 ,  60) (RE(L) ,  L  -  1,  NR) 

60  FORMATC  INTERIOR  FIELD  AMPLITUDE’ /(IX,  4EI4.7)) 
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The  above  write  statements  are  labeled  the  first,  second,  third,  .... 
eighth  write  statements.  The  data  of  the  first  write  statement  are 
merely  input  data.  The  data  of  the  second  and  third  write  statements 
are  defined  by 

(BJ(NP),  NP  -  JJ,  JJN)  *  (J  (ka),  n-0,1 .  Nl-1)  (17) 

n 

(BY (NP) ,  NP  -  JJ,  JJN)  =  (Y  (ka),  n  -  0,  1 . Nl-1)  (18) 

n 

where  J  (ka)  and  Y  (ka)  are,  respectively,  the  Bessel  functions  of  the 
n  n 

first  and  second  kinds  of  order  n  and  argument  ka.  The  data  of  the  fourth 
and  fifth  write  statements  are  defined  by 


(SE(J) ,  J-l,  NN)  -  (|v“|  ,  j-1,2 . NN)  0 

(S0(J) ,  J-l,  NN)  -  (|V°|  ,  j-1,2 . NN)  (5 

6  O 

where  appears  in  [1,  eq.  (7)J  and  V\  appears  in  [1,  eq.  (24)].  If 

6  O 

NN  <  4,  then  and  are  truncated  at  j  *  NN  because  they  are  given  in 

Ag  AQ 

terms  of  and  by  [1,  eqs.  (53)  and  (54)]  and,  as  ind:  ated  earlier, 

Ag  A.Q 

V.  and  V,  are  truncated  at  1  =  NN. 

j  j 

In  the  sixth  write  statement,  U2  is  the  complex  number  whose  real 
and  imaginary  parts  are  those  of  E^(-M)  at  the  center  of  the  cylinder. 


and  SI  is 


|E^(-M)  |  at  the  center  of  the  cylinder.  This  E^(-M)  is 


calculated  from  [1,  eq.  (68)].  In  the  seventh  write  statement,  BJ(L)  is 

the  magnitude  of  the  aperture  field  E  at  p  «  a  and  <{>-<]>  where 

Z  L 

2(L  -  1)  4> 


In  the  eighth  write  statement,  RE(L)  is  |E  (-M) |  at  ko  *  ko  along  the 

Z  L 

line  that  runs  from  the  center  of  the  cylinder  to  the  center  of  the  aper¬ 
ture.  Here, 


where 


Y*  *  l  E .  (EF) 
ij  n-0  Jn  in 


jj=2,3,< 


H(2)(ka) 

n 


e  e  e 

s,  +  c  s, 

_in _ j In 

J  (ka) 
n 


e  e  e 

S.  +  C  4  S- 

Jn  J  -  In 

J  (ka)  Y  (ka) 
n  n 


0  ,  j-1 


n-0,1,. • • ,N1-1 


j  F  ,  n  *  Nl,  Nl+1 , . . . ,N2 


,  n*0,l , . . . ,p-l ,p+l . Nl-1 


,  n=Nl,  Nl+1 . N2 


,  J-2,3,4 


Assuming  that  (23)  holds,  we  approximate  the  matrix  elements 
[1,  eqs.  (47),  (48),  and  (50)]  by 


F  S°  Nl-1  F°  S°  N2  F°  S, 

_j£ _ iE_  +  j0(ka)  Y  _ ^  r  m  In 

(2)  VK  i  T  ty,  ^  u(2)/.  \  in  J  (ka)  Y  (ka) 

H  (ka)  r  n=l  J  (ka)  H  (ka)  n=Nl  n  n 

p  Lfp  n 


i-1,2,3,4 


N1_1  F°  (S°  +C°S°  )  N2  F°  (S°  +C°S°  ) 

y  in  jn  J  In  +  .  y  in  jn  j  In 

L,  t  \  J  J  (ka)  Y  (ka)  ’ 

n=l  J  (ka)H  (ka)  n=Nl  n  n 

,  n  n 

n#p 


i-1,2,3,4 


j-2,3,4 


i*l »2 ,3,4 


Nl-1  Fj  sin  (na) 

2  y  _ _ in - 

z  4  (2) 

n-1  Hu;(ka) 
n 


where 


J° (ka) 
P 


1  .  p  -  0 


J  (ka)  ,  p  >  0 
P 


and  [1*  eq.  (35)] 


0  ,  p  -  0 


■f  ’  P>° 

lp 

Expressions  (33)  and  (34)  are  recasC  as 


N2 


+  J>»>  I  Eu<0F)ln  '  “1>2’3>4 


il  lp'  ip 


n«0 

nj*p 


N2 

,o  r  _o 


■ii-J0Von*»  1 

n^p 


where 


(OF) 


Jjn 


23 


0  ,  P  -  0,  j-1,2,3,4 


0  ,  p  >  o,  j«l 


A 


,  P  >  0,  j-2,3,4 


The  terms  for  which  n-0  do  not  contribute  to  (38)  and  (39)  because 

S°  -  F°  *  0  [1,  eqs.  (27)  and  (49)].  The  terms  for  which  n-0  have  been 
j  o  io  * 

added  to  make  (38)  and  (39)  similar  to  (28)  and  (29)  and  thus  easier  to 


program. 


For  convenience,  [1,  eqs.  (40),  (46),  and  (53)]  are  repeated 


e*e  -*e 
Y  V  -  I 


O^O 

Y  V  -  I 


V*-Jp(ka)V*+  I  j 

1-2  V  (45 

vj  "  ^  ‘  j"2’3’4  J 

With  J°(ka)  given  by  (36)  and  C°  by  (42),  [1,  eqs.  (54)  and  (55)]  combine: 


V°  -  Jp(ka)V°  +  l  C°V° 


V°  -  V°  ,  j-2,3,4 


Expression  [1,  eq.  (68)] for  the  field  E^(-M)  at  the  center  of  the 


cylinder  is 


6  ^  6 
S  V 
10  1 


j, 


p  -  o 


,  p  +  o 
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The  aperture  field  E  is  the  cj)  component  of  M  of  [1,  eq.  (4)]: 

z  — 

E  -  Me  +  M°  (48) 

z 

Substituting  [1,  eqs.  (7)  and  (24)]  into  (48),  using  [1,  eqs.  (8)  and 
(25)],  and  evaluating  E  at  4>  ■  {<L  ,  L  »  1,2 . NA}  where  <L  is  given  by 

Z  L  L 


(21),  we  obtain 


4  <j>. 


[Ez]L  “  [  (VQ  vj  +  (^)2J"2]  / 1  '  (^)2  ’  L“1*2*-‘-*NA  (49> 

Along  the  line  from  the  center  of  the  cylinder  to  the  center  of  the 
aperture,  4)  ■  0  so  that  expression  [1,  eq.  (62)]  for  the  interior  field 


E  (-M)  reduces  to 
z  — 


,  °°  B  J  (ko) 

Ez("^  ”  l  j  (klT 

Z  n-0  nUa} 


where  is  given  by  [1,  eqs.  (63)  or  (65)]. 


VJ  ' "  * 


jyka)  V*S®p  ,  n-p 

Truncating  the  summation  in  (50)  at  n  -  N2,  substituting  (51)  into  (50), 
and  then  letting  kp  *  {kp^.  L-l ,2 , . . . ,NR)  where  kp^  is  given  by  (22), 
we  have 


4  N2  S,  J  (kPT) 

'.e  .  r  „e  ,  r  jn  n  L  , 


[Ez(-«)]L  "  SlpVkPlA  +  jI1  Vj  (JmQ  VCka)  >  *  (52) 

n^p 


Expression  (52)  is  recast  as 


eVm)].  -  (RP).tf®  +  l  (RN)  V®  ,  L  -  1,: 
z  -  L  L  1  j-1  jL  J 


where 
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(RP)L  “  Slp  Jp(kpL> 


Nl-1  S®  J  (kp.)  N2  S®  J  (kp.) 
r  jn  n  L  ,  r  inn  L 

l  T  (V  +  I  T  Hca} 


J  (ka) 
n 


J  (ka) 
n 


Expression  (55)  was  obtained  by  assuming  that  (23)  holds.  We  plan  to 
use  Debye's  asymptotic  expansions  for  the  Bessel  functions  in  the  second 
summation  in  (55). 

If  L-l,  then  kp^  =  0  and  (53)  reduces  to 


[Eb(-M)] 
z  —  i 


P  AP 

S  V 
10  1 


4  S. 


t  P  *  0 


^  vj  •  >* 0 


Therefore,  any  difference  between  the  computed  values  of  the  field  (47)  at 
the  center  of  the  cylinder  and  the  field  (53)  at  L  *  1  is  due  to  roundoff 


error. 


If  L  ■  NR,  then  kp^  -  ka  and  (52)  reduces  to 


*  SlPVka)?!  +  £  v]  ( jQ  S]n> 

n^p 


Using  (45)  to  express  J^(ka)V^  in  terms  of  {V^,  j=l,2,3,4},  we  obtain 

J  (ka)vj  -  V®  -  l  C®V® 

P  1=2 

In  view  of  (32),  substitution  of  (58)  into  (57)  gives 

b  *  e  Nv2  SG  ) 

CEz(^)]NR  *  l  V1(  l  V 
2  NR  j-1  j  n-0 

Equating  [1,  eqs.  (8)  and  (12)]  at  <t  ■  0,  we  obtain 


n*0  J 


1  ,  j-1 


0  ,  j-2,3,4 


m 

a$ 


bvy 


I .v*v 


V 


w.-.vs-v 


26  . 


(61) 


(62) 


If  NA  Is  odd  so  that  *  0  when  L  ■  (NA+l)/2,  then  (49)  gives 

[El.  -  V*  ,  L  -  (NA+l)/2  (63) 

Z  L  1 

Noting  that  (62)  and  (63)  have  the  same  right-hand  side  and  recalling 
that  (61)  was  used  to  obtain  (62),  we  deduce  that  any  difference  between 
the  computed  values  of  the  aperture  field  (49)  at  the  center  of  the  aper¬ 
ture  and  the  interior  field  (53)  at  kp^  ■  ka  is  due  to  the  error  in  truncating 
the  summation  with  respect  to  n  in  (55),  the  error  in  calculating  {S^},  and, 
of  course,  roundoff  error. 

Having  presented  relevant  formulas,  we  begin  to  verbally  describe  the 

0 

main  program.  In  the  equations  that  were  presented  in  this  section,  Y  and 

Y°  are  4x4  matrices,  v  ,  I6,  $°,  and  T3  are  4x1  column  vectors,  and 

{Vj,  V°,  j“l,2,3,4}  are  defined  by  (45)  and  (46).  The  equations  that  were 

programmed  are  those  that  were  presented  in  this  section  with  the  order  4 

e  o 

replaced  by  NN,  that  is,  Y  and  Y  are  replaced  by  the  NN x NN  submatrices 

^6  0  0  6  O 

in  their  upper  left-hand  corners,  and  Vj ,  Ij ,  V  ^  ,  1^ ,  ,  and  are  re¬ 

tained  only  for  j  £  NN.  We  describe  the  main  program  by  relating  variables 
therein  to  quantities  that  appear  in  the  previous  text  of  this  section.  What 
follows  is  clarified  in  Table  1. 

Line  27  and  DO  loop  61  put 


a 


in  BJ(n+l)  for  n-0,1, . . . ,N1-1.  Equation  (64)  can  be  obtained  from  [3, 


formula  9.1.10].  In  DO  loop  18,  line  38  puts  J  (kp  )  in  BJ (n+l+(L-l)*Nl) 

n  Li 


for  n-0,1, ... ,N1-1  where  kp.  is  given  by  (22).  Line  38  also  puts  Y  (kpT ) 

L  n  L 


in  BY (n+l+(L-l)*Nl)  for  n-0,1 . Nl-1  although  only  {Y  (kp  1  -  Y  (ka), 

n  nk  n 


n-0,1,  ..., Nl-1 }  is  needed. 


Table  1.  Key  lines  in  the  main  program,  the  variables  read  in,  defined, 
or  incremented  therein,  the  corresponding  quantities  in  the 
text,  and  the  numbers  of  the  equations  where  these  quantities 
appear . 

Line  Variable(s)  in  the  main  .^Quantity  in  the  text  |  Equation 

Number  program 


Number 


|  Equation 
I  number (s) 


IP  and  X 

P  and  ALP 

BJ(1) 

BJ(NP) 

XR(1) 

XR(L) 


Ip  and  ka 


and  a  in  radians 


J  (0) 

o 


W°> 


f  BJ(n+JJ)  ,  n-0,  1 . Nl-1)  {.I  (k0,),  n-0,1 . Nl-1)  26,30,31,54,55 

n  L 

(BY (n+JJ)  ,  n-0,  1.....N1-1}  (Y  (kpT  )  ,  n-0,1 . Nl-l}  ',26,30 

n  L  i 


{ SE(J) ,  J-l ,2,3,4} 
{S0(J),  J-l, 2, 3, 4} 
BJPE 
CE(1) 

CE(J) 


Nl-1} 

{ J  (kp, 

n 

Nl-1} 

(Yn(kp, 

(s°  . 

|  JP 

I 

1 J  (ka) 

-j 


Variable (s)  in  the  main 
Number  I  program 


Quantity  in  the  Text 


125 

{BJJ(L),  L-1,2,...,NR} 

{J  (kp  ) /J  (ka)} 
n  L  n 

128.131 

E 

EJ 

129,132 

0 

0 

J 

135 

YE(IY) 

nth  term 

136 

Y0(IY  3 

nth  term 

143 

RN(JR) 

nth  term 

148,152 

{VE(J  )  ,  1-1,2 . NN} 

{V®,  1-1,2,. ...NN} 

149,154 

{VO(I),  1-1,2 . NN} 

{V°,  1-1,2 . NN} 

163 

V0(1) 

o  «— < 
> 

164 

U3 

165 

VE(1) 

+1 

167 

SE(  J) 

lvj ! 

168 

S0(J) 

|vj ! 

173 

SE(J ) 

{sj<r 

175 

U2 

se  ve 

S10  1 

181 

U2 

188 

XX 

V*o 

192 

U2 

195 

BJ  (L) 

K'iJ 

201 

U2 

^  P 

<»Vl 

204 

U2 

206 

RE(L) 

1  [  E^  ( -M) 1 L | 

Equation 
number (s) 


28  or  29 
38  or  39 


3' 


V 


Line  61  puts  S®  of  (32)  in  SE(J)  for  J-1,2, 3,  and  4.  Line  61  also 
Jp 

puts  S°  of  (42)  in  SO(J)  for  J-1,2, 3,  and  4.  Lines  62  to  67  put  J  (lea)  of 

Jp  p 

(28)  in  BJPE  and  {C®,  J-1,2 . NN>  of  (32)  in  (CE(J),  J-1,2,. . . ,NN}.  Lines 

68  to  79  put  J°(ka)  of  (33)  in  BJPO  and  (C°,  J-1,2, ... ,NN}  of  (42)  In 
(CO(J),  J-1,2 . NN}. 


Being  sumnations  with  respect  to  n,  the  matrix  elements  I®  of  (26), 
of  (28)  and  (29),  1°  of  (35),  Y°j  of  (38)  and  (39),  and  (RN)^  of  (55) 
are  accumulated  in  CUE(I),  YE(IY) ,  CUO(I),  YO(IY),  and  RN(JR) ,  respectively, 
where 


IY  -  I  +  (J-1)*NN 

» 

JR  -  J  +  (L-1)*NN 

In  DO  loop  15, 


(65) 


N  -  NP-1 


(66) 


where  NP  is  the  index  of  DO  loop  15.  The  terms  for  which  n-N  are  added  to 

CUE(I)  and  CUO(I)  in  DO  loop  15.  If  N  j*  p,  then  the  terms  for  which  n-N  are 

added  to  YE(IY),  YO(IY),  and  RN(JR)  in  DO  loop  15.  If  N-p,  then  S®  (EF)_  is 

lp  Ip 

added  to  YE(I),  S?  (0F)T  is  added  to  YO(I),  and  (RP)T  of  (54)  is  put  in  RP(L) 
lp  lp  L 

in  DO  loop  15.  More  detail  on  the  workings  of  DO  loop  15  is  presented  in 
the  following  four  paragraphs. 

Line  85  puts  S®  of  (31)  in  SE(J),  and  S°  of  (41)  in  SO(J)  for  J-1,2, 3, 
jn  Jn 

and  4.  Line  85  also  puts  F®n  of  (26)  and  (30)  in  FE(I),  and  F°n  of  (35)  and 
(40)  in  F0(I)  for  1-1,2, 3  and  4. 

If  N<N1,  then  lines  87  to  118  are  executed.  Line  88  puts  J  (ka) 

n 

in  BJN,  line  91  puts  1/H^(ka)  in  U2,  lines  93  and  94  put  c  cos  (na)  in 

n  n 

CSA,  and  line  95  puts  2sin(na)  in  SNA.  In  DO  loop  40,  (EF)^^  of  (30)  is 


put  in  EF(I),  (OF)^  of  (40)  is  put  in  0F(I),  the  nth  terra  of  (26)  is 

added  to  CUE(I),  and  the  nth  term  of  (35)  is  added  to  CUO(I).  If  N»p, 

then  DO  loops  41  and  53  are  executed.  In  DO  loop  41,  sf  (EF),  of  (28) 

Ip  Ip 

is  added  to  YE(I),  and  S°  (OF)  of  (38)  is  added  to  Y0(I).  In  DO  loop  53, 

lp  Ip 

(RP)^  of  (54)  is  put  in  RP(L).  If  N  p,  then  DO  loop  55  is  executed.  In 

DO  loop  55,  the  factor  J  (koT)/J  (ka)  in  the  first  summation  in  (55)  is  put 

n  L  n 

in  BJJ(L). 

If  N  N1 ,  then  lines  120  to  125  are  executed.  Line  120  puts 

(ka)  Yn(ka)  of  (31)  and  (41)  in  BJN.  DO  loop  51  puts  jF®n  of  (30)  in 

EF(I)  and  jF°  of  (40)  in  0F(I).  In  line  125,  the  factor  J  (kpT)/J  (ka) 

In  n  L  n 

in  the  second  summation  in  (55)  is  put  in  BJJ(L)  for  L=*l  ,2 , . .  .  ,NR. 

If  N  i  p,  then  lines  126  to  145  are  executed.  Immediately  before 

these  lines  are  executed  (EF)_  of  (30)  will  reside  in  EF(I),  and  (OF),. 

In  In 

of  (40)  will  reside  in  0F(I).  Upon  entry  into  DO  loop  44,  the  quantity 
Ej  defined  by 


E. 

J 


resides  in  E  and  the  quantity  0 


J 


J  (ka)  Ee  , 
P  jn 


defined  by 


j  “  1 
j  *  1 


(67 


J°(ka)  E?  ,  J  -  1 
P  jn’ 


(68 


»  j  l4  1 


resides  in  0.  Here,  E^  is  given  by  (31)  and  E°^  by  (41).  In  DO  loop  44, 


the  nth  term  of  the  summation  in  (28)  or  (29)  is  added  to  YE(IY),  and 
the  nth  term  of  the  summation  in  (38)  or  (39)  is  added  to  YO(IY).  In 


nested  DO  loops  56  and  57,  the  nth  term  of  the  summation  in  (RN)^  of 


(55)  is  added  to  RN(JR). 


Lines  147  to  154  put  the  elements  of  the  solution  V  of  (43)  in 
VE  and  the  elements  of  the  solution  V°  of  (44)  in  VO.  In  DO  loop  45,  V^ 
of  (45)  is  accumulated  in  U2  and  V°  of  (46)  is  accumulated  in  U3.  The  series 
in  (45)  and  (46)  were,  as  Indicated  in  the  paragraph  preceding  (64),  summed 
for  1  up  to  NN  rather  than  4.  After  lines  163-165  have  been  executed,  the 

-fg  Ag 

elements  of  V  will  reside  in  VE,  the  elements  of  V  will  reside  in  VO,  and  V^ 
will  be  in  U3.  In  DO  loop  47,  |v®|  is  put  in  SE(J) ,  and  |V°|  is  put  in  SO(J). 

The  rest  of  the  main  program,  lines  173  to  210,  serves  to  calculate 

the  field  (47)  at  the  center  of  the  cylinder,  the  aperture  field  (49),  and 
the  field  (53)  along  the  line  from  the  center  of  the  cylinder  to  the  center 
of  the  aperture.  The  series  in  (47),  (49),  and  (53)  are  summed  with  i  running 
from  1  to  NN  rather  than  from  1  to  4. 

Line  173  puts  S*Q  in  SE(J)  for  J-1,2,3,4.  If  p-0,  then  line  175 
puts  of  (47)  in  U2.  If  p  +  0,  then  DO  loop  48  is  executed  and 

line  181  puts  [Eb(-M)]p-0  of  (47)  in  U2. 

The  index  L  of  DO  loop  62  obtains  the  subscript  L  in  (49).  Line  188 
puts  (<tT/4>  )  in  XX.  DO  loop  63  accumulates  [E  ]T  of  (49)  in  U2.  Line  195 

puts  |[E  L  j  of  (49)  in  BJ(L). 

The  index  L  of  DO  loop  58  obtains  the  subscript  L  in  (53).  Line  201 
puts  (RP)^V®  of  (53)  in  U2.  DO  loop  59  adds  to  U2  the  summation  with  respect 
to  j  in  (53).  Line  206  puts  |(Eb(-M)l  |  of  (53)  in  RE(L) . 

Sample  input  and  output  data  appear  after  the  listing  of  the  main 
program.  The  sample  input  data  are  for  a  perfectly  conducting  circular 
cylindrical  shell  of  electrical  radius  ka  •  2  with  a  slot  of  half  angular 
width  4>q  *  1.25°  ■  0.02181662  radians  centered  at  4)  ■  0°  illuminated  by  a 
TM  plane  wave  that  comes  from  the  direction  in  which  ■  0°.  The  last 


number  in  the  sixth  from  the  last  row  of  the  sample  output  data  is  exactly 
the  same  as  the  first  number  in  the  last  row.  This  is  expected  because  ex¬ 
pression  (56)  is  identical  to  expression  (47)  when  p-0. 

We  ran  the  main  program  with  IP  of  the  input  data  changed  from 
0  to  1.  The  first  four  significant  figures  of  each  number  in  the  last  six 
lines  of  the  resulting  output  data  were  the  same  as  those  that  were  obtained 
with  IP  -  0.  According  to  (56)  and  (47),  the  last  number  in  the  sixth  from  the 
last  line  of  the  output  data  may,  due  to  roundoff  error,  be  slightly  dif¬ 
ferent  from  the  first  number  in  the  last  line  if  IP^O.  These  two  numbers 
were  identical  to  each  other  in  the  output  obtained  with  IP*1. 

Note  that  the  third  number  in  the  fourth  from  the  last  line  of  the 
sample  output  data  is  slightly  different  from  the  last  number  in  the  last 
row.  Reasons  for  this  are  given  in  the  paragraph  that  contains  (63). 
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1C  LISTING  OF  THE  RAIN  PROGRAM 

2 C  INPUT  DATA  ARE  READ  FROH  THE  FILE  HAUTZ3.DAT  ON  UNIT  20. 

3 C  OUTPUT  DATA  ARE  WRITTEN  ON  THE  FILE  HAUTZ6.DAT  ON  UNIT  21. 

4 C  THE  SUBROUTINES  BES.  BESJ,  BESJY.  BESJJ.  SFEO. 

5  C  DECO HP .  AND  SOLVE  ARE  USED. 

6  COMPLEX  YE( 16) • Y0( 16) . CUE(4) • CUO (4) • U3. U. U2. EF (4) . OF (4) 

7  COMPLEX  VE(4).V0(4) 

8  DIMENSION  Y0(36) ,BJ (220) .BY (220) . XR( 1 1 ) . RP ( 1 1 ) . RN(44) . RE( 1 1 ) 

9  DIMENSION  BJJ(11).SE(4),S0(4),FE(4),F0(4).CE(4),C0(4). IPS<4> 

10  OPEN (UN IT-20,  FILE- 'HAUTZ3.DAT ’.STATUS- 'OLD* ) 

11  OPEN (UN IT-21. FILE- 'HAUTZ6.DAT' . STATUS-’OLD’ ) 

12  READ(20.25)  (YOU).  1-1.36) 

13  25  FORMAT (5E14. 7) 

14  WRITE(21. 26)  (YOU),  1-1.38) 

15  26  FORMAT ( '  TO'/ (1X.5E14.7) ) 

16  READ (20. 28)  N1.N2.N3.NA.NR 

17  28  FORHAT( 13. 15.313) 

16  VRITE(21.29)  N1.N2.N3.NA.NR 

19  29  FORMAT (’  Ml-'. 13.’.  M2-'. 15.’.  N3-*.I3.».  NA-'.I3.’.  NR»».13) 

20  M1H-N1-1 

21  DO  31  JV-1.N3 

22  READ (20. 32)  NN.  IP. IB. X. P. ALP 

23  32  F0RHAT(3I3.3E14.7) 

24  WRITE(21 • 27)  NN. IP. IB.X.P. ALP 

25  27  FORMAT ('  NN«’.I3.*.  IP-', 13.'.  !B«'.I3/ 

26  1*  X— ' , E14. 7, ' ,  P-' , E14. 7, * .  ALP-’.E14.7) 

27  BJ ( 1)«1. 

28  00  61  NP-2.N1 

29  BJ (NP)-O. 

30  61  CONTINUE 

31  XR(1)*0. 

32  JJ-1 

33  FNR-NR-1 

34  00  16  L-2.NR 

35  XR(L)«(L-l)/FNR»X 

36  IF(L.EQ.NR)  XR(L)-X 

37  JJ-JJ+N1 

38  CALL  BES( Y0.N1M, XR(L) . BJ ( JJ ) . BY( JJ ) ) 

39  18  CONTINUE 

40  JJN-JJ+N1-1 

41  WRITE(21.33) (BJ(NP).NP-JJ. JJN) 

42  33  FORMAT ( ’  BJ'/ (1X.5E14.7) > 

43  WRITE(2l , 19) (BY(NP) .NP-JJ. JJN) 

44  19  FORMAT ( »  BY' / (IX, 5E14. 7) ) 

45  NN2«NN«NN 

46  DO  35  I-1.NN2 

47  YE( I )»0. 

48  Y0( I )-0. 

49  35  CONTINUE 

50  DO  36  1*1. NN 

51  CUEU)«0. 

52  CU0(I)-O. 

53  36  CONTINUE 


54 

NRN-NR*NN 

55 

00  52  J-l.NRN 

56 

RN<J>«0. 

57 

52 

CONTINUE 

58 

N1R*N1«(NR-1) 

59 

IPP-IP+1+N1R 

60 

I F < IB.EQ.O)  BJ ( IPP) *0. 

61 

CALL  SFEQ( IP.P.SE.SO.FE.FO) 

62 

BJPE-BJtIPP) 

63 

CEC 1  )-0. 

64 

IFtNN.EQ. 1)  GO  TO  34 

65 

DO  37  J-2.HN 

66 

CE ( J ) «-SE ( J ) /SE  < 1 ) 

67 

37 

CONTINUE 

68 

34 

IF( IP.NE.O)  GO  TO  21 

69 

BJPO-1. 

70 

DO  38  J-l.NN 

71 

C0<J)»0. 

72 

38 

CONTINUE 

73 

GO  TO  23 

74 

21 

bjpo-bjpe 

75 

C0<1>»0. 

76 

IFtNN.EQ. 1)  GO  TO  23 

77 

00  22  J-2.NN 

78 

COtJ)— S0(J>/S0<1> 

79 

22 

CONTINUE 

80 

23 

U3-U..0.  ) 

81 

U-<0. , 1. ) 

82 

N2P-N2+1 

83 

DO  15  NP-1 . N2P 

84 

N-NP-1 

85 

CALL  SFEO(N.P.SE.SO.FE.FO) 

86 

IF(N.GE.Nl)  GO  TO  50 

87 

NPN-NP+N1R 

88 

BJN-BJ(NPN) 

89 

Sl-BJN 

90 

IF(ABStSl) . LE. l.E-10)  Sl-0 

91 

U2- 1 . / CHPLX ( S 1 . -B Y  <  NPN  > ) 

92 

ALN-ALP»N 

93 

CSA-2. oCOS(ALN) 

94 

IF(N.EQ.O)  CSA-1. 

95 

SNA-2. *SIN<ALN> 

96 

DO  40  I-l.NN 

97 

EF< I )»FE< I ) *U2 

98 

OF ( I )-F0 ( I ) *U2 

99 

CUE ( I ) »CSA«EF ( I ) »U3*CUE  < I ) 

100 

CUO  < I ) «SNA*OF  t I ) »U3*CU0 < I > 

101 

40  CONTINUE 

102 

U3-U3»U 

103 

IFtN.NE. IP)  GO  TO  20 

104 

DO  41  I-l.NN 

105 

TE ( I ) -SE ( 1 ) *EF ( 1 ) ♦YE ( I > 

106 

YO  < I ) -SO  < 1 ) *0F  < I > ♦YO ( I ) 

.rfV  ■' V' -*»  -r.  V  v. 


107 

41 

CONTINUE 

108 

1J-NP 

109 

00  53  L«1.NR 

no 

RP(L)*SE( 1 ) »BJ ( I J) 

111 

I J*l J+Nl 

112 

53 

CONTINUE 

113 

GO  TO  15 

114 

20 

JJ»NP 

115 

00  55  L-l.NR 

116 

BJJ  <L)«BJ ( JJ)/BJN 

117 

JJ-JJ+N1 

118 

55 

CONTINUE 

119 

GO  TO  54 

120 

50 

CALL  BESJY(N.X.BJN) 

121 

DO  51  I»1.NN 

122 

EF( I )*FE( I )»U 

123 

OF ( I >*F0 ( I ) »U 

124 

51 

CONTINUE 

125 

CALL  BESJJ (NfNRtXRf BJJ) 

126 

54 

IT=»0 

127 

DO  42  J*1 • NN 

128 

E*(SE( J ) +CE(  J ) »SE( 1 ) ) /BJN 

129 

0-(S0(J)+C0<J)*S0U>  >/BJN 

130 

IF(J.NE.l)  GO  TO  43 

131 

E-BJPE«E 

132 

0*BJP0»0 

133 

43 

DO  44  I-l.NN 

134 

ir»ir+i 

135 

• 

YE(IY)-E»EF(I)+TE(IY) 

136 

Y0( IY)*0»0F( I )+Y0( IY) 

137 

44 

CONTINUE 

138 

42 

CONTINUE 

139 

JR*0 

140 

DO  56  L*  1  •  NR 

141 

00  57  J«1.MN 

142 

JR-JR+1 

143 

RN ( JR) *BJ J (L) *SE ( J ) +RN  <  JR) 

144 

57 

CONTINUE 

145 

56 

CONTINUE 

146 

15 

CONTINUE 

147 

IFINN.NE. 1 )  GO  TO  17 

148 

VE ( 1 ) "CUE ( 1 ) / YE ( 1 ) 

149 

VO  ( 1 ) *CU0 ( 1 ) / YO ( 1 ) 

150 

GO  TO  16 

151 

17 

CALL  DEC0I1P  ( NN  •  I  PS  >  YE  > 

152 

CALL  SOLVEINN. IPS. YE.CUE. VE) 

153 

CALL  DECOHPINN, IPS.YO) 

154 

CALL  SOLVEINN. IPS.YO. CUO. VO) 

155 

16 

CE(1)«BJPE 

160  U2»CE(I)*VE<I)+U2 

161  U3-C0(I)*V0(I)+U3 

162  45  CONTINUE 

163  V0(1)«U3 

164  U3*VE  < 1 ) 

165  VE(1)-U2 

166  DO  47  JM.NN 

167  SE<J>«CABS<VE(J>> 

168  SO(J)-CABS(VO(J)> 

169  47  CONTINUE 

170  URITE<21.24> <SE< J) . J-l.NN) 

171  VR1TE(21.24) (SO(J).J-l.NN) 

172  24  F0RHAT(1X*4E14.7) 

173  CALL  SFEO  ( 0 .  P .  SE .  SO .  FE .  FO ) 

174  IF( IP.NE.O)  GO  TO  46 

175  U2«SE(l)»U3 

176  GO  TO  39 

177  46  U2«0. 

176  DO  46  J*1  •  NN 

179  U2«SE ( J ) * VE ( J ) +U2 

180  46  CONTINUE 

181  U2“U2/BJ (N1R+1 ) 

182  39  S1*CABS(U2) 

183  «RITE(21.49)  U2.S1 

184  49  FORHAT ( *  £■' • 2E14.7. * .  ABS<E)>* > E14. 7) 

185  FNA*2./ <NA“1) 

186  DO  62  L«1.NA 

187  U2-0. 

188  XX»-1.+<L-1)*FNA 

189  X2-XX»XX 

190  S1=»SQRT<1.-X2> 

191  00  63  J«1  •  NN 

192  U2-S1*(XX*V0(J)«VE(J) )«U2 

193  S1-S1»X2 

194  63  CONTINUE 

195  BJ(L)*CABS(U2) 

196  62  CONTINUE 

197  URITE(21.64) (BJ(L).L-l.NA) 

198  64  FORHAT ( *  APERTURE  FIELD  AHPLITUDE' / < IX. 4E14. 7) > 

199  JR*0 

200  DO  56  L-l.NR 

201  U2«RP(L)»U3 

202  DO  59  J*1.NN 

203  JR«JRM 

204  U2*RN  < JR) »VE  <  J ) +U2 

205  59  CONTINUE 

206  RE(L)*CABS(U2) 

207  58  CONTINUE 

206  WRITE(21.60)<RE(L).L»1.NR) 

209  60  FORHAT ( *  INTERIOR  FIELD  AHPLITUDE*/ ( 1X.4E14. 7) ) 

210  31  CONTINUE 

211  STOP 

212  EHD 


r>  o 
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C  LI ST IMG  OF  THE  IMPUT  DATA  FILE  HA0T23.DAT 

-0. 3072582E+04  0.7368758E+04-0.6085100E+03  0. 1710234E+02-0. 2271001E+00 
0. 16001 71 E-02-0. 5961 089E-05  0.9545773E-08  0. 4163150E+05  0. 342021 1E+03 
0. 1 OOOOOOE+O 1 -0 . 6024727E+04  0. 1613512E+04-0. 7532210E+02  0. 1402590E+01 
-0 . 1 275602E-0 1  0 . 5832787E-04-0 . 1 1 07698E-06  0 . 3072946E+05  0 . 2886431 E+03 
0. 1 OOOOOOE+O 1  0.9999999E+00-0. 1097659E-02  0. 2461455E-04  0. 1 OOOOOOE+O 1 
0 . 1 829893E-02-0 .3191 328E-04-0 . 1 562498E-0 1  0.1 4270 79E-03-0 . 5937434E-05 
0 . 4687498E-0 1-0.1 998720E-03  0. 7317495E-05  0. 6366198E+00  0. 7853982E+00 
0.2356194E+01 
2010000  1  5  3 

4  0  1  0 . 2000000E+0 1  0. 2181662E-01  0. OOOOOOOE+OO 

LISTING  OF  THE  OUTPUT  DATA  FILE  HAUTZ6.DAT 
TO 

-0. 3072582E+04  0. 7368758E+04-0. 6085100E+03  0. 1710234E+02-0. 227100 1E+00 
0. 16001 7 1 E-02-0. 596 1089E-05  0. 9545773E-08  0. 4163150E+05  0. 342021 1E+03 
0. 1 OOOOOOE+O 1-0. 6024727E+04  0. 161 35 12E+04-0. 75322 10E+02  0. 1402590E+01 
-0. 1275602E-01  0. 5832787E-04-0. 1 107698E-06  0. 3072946E+05  0.2886431E+03 
0. 1 OOOOOOE+O 1  0. 9999999E+00-0. 1097659E-02  0. 2461455E-04  0. 1 OOOOOOE+O 1 
0. 1 829893E-02-0 .3191 328E-04-0 . 1562498E-01  0. 1427079E-03-0.5937434E-05 
0 . 4687498E-0 1-0.1 998720E-03  0 . 73 1 7495E-05  0 . 6366 1 98E+00  0 . 7853982E+00 
0.2356194E+01 

Nl-  20.  N2-10000.  N3»  1.  HA*  5.  NR*  3 

NN«  4.  IP*  0.  IB*  1 

X*  0. 2000000E+01 •  P-  0. 2181662E-01 .  ALP*  0. OOOOOOOE+OO 
BJ 

0. 2238908E+00  0.5767248E+00  0. 3528340E+00  0. 1289432E+00  0. 3399571E-01 
0. 7039629E-02  0. 1202429E-02  0. 174944 IE-03  0. 2217955E-04  0. 2492343E-05 
0. 2515386E-06  0. 2304285E-07  0. 1932695E-08  0. 1494942E-09  0. 1072946E-10 
0. 7183016E-12  0.4506005E-13  0. 2659308E-14  0. 1481737E-15  0. 7819242E-17 
BY 

0.51 03757E+00-0. 1 070324E+00-0 . 6 1 7408 1 E+00-0. 1 1 27784E+0 1 -0 . 2765943E+0 1 
-0 . 9935989E+0 1 -0 . 469 1 40 1 E+02-0 .2715481 E+03-0 . 1 853922E+04-0 . 1 455983E+05 
-0. 1291846E+06-0. 1277286E+07-0. 1392096E+08-0. 1557742E+09-0. 2141 144E+10 
-0 . 298 1 024E+ 1 1 -0 . 4450 1 25E + 1 2 -0 . 7090389E ♦ 1 3-0 . 12009 1 6E ♦ 1 5-0 . 2 1 54558E ♦ 1 6 
0. 4655501E-01  0. 1714323E-01  0. 6358480E-01  0. 5726649E-01 
0. OOOOOOOE+OO  0. OOOOOOOE+OO  0. OOOOOOOE+OO  0. OOOOOOOE+OO 
E*-0 .11211 77E-02-0 . 272802 1 E-03 .  ABS ( E ) *  0.11 53889E-02 
APERTURE  FIELD  AHPLITUDE 

0. OOOOOOOE+OO  0. 4136272E-01  0. 4655501E-01  0. 4136272E-01 
0. OOOOOOOE+OO 
INTERIOR  FIELD  AHPLITUDE 
0. 1 153889E-02  0. 1597396E-02  0. 4656652E-01 
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